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ABSTRACT 


In  this  work  the  prohlem  of  a  supercavitating  hydrofoil  in 
linear  theory  is  stated  in  terms  of  integral  equations  of  unknown 
vortex  and  cavity  source  distributions.  The  general  problem  is 
decomposed  into  the  camber,  the  thickness  and  the  angle  of  attack 
problems.  For  each  problem  the  solution  is  expressed  in  terms  of 
integrals  of  the  slope  of  the  lower  surface  of  the  hydrofoil  weighted 
by  known  functions.  The  numerical  scheme  to  compute  those  integrals 
has  been  proved  very  accurate  and  insensitive  to  the  variables  of  the 
problem.  As  an  application,  tables  and  some  cavity  plots  are  given  for 
the  NACA  a  ■  0.8  meanline  series  and  the  NACA  Four-digit  Wing  Sections 
thickness  forms. 
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NOMENCLATURE 


C 

C('x) 

C 


/ 

0) 

F(t) 

m 

& 

^Cx) 

IM 

m 

3h) 

K 

Kfz) 


Constant  defined  as  Ct  s 
Constant  defined  as  b 
Integral  defined  in  (8.4) 

Chord  of  the  hydrofoil 

Camber  distribution  of  the  foil  and  cavity  (see  Figure  1) 

Constant  defined  in  (9.5) 

Drag  coefficient 

Lift  coefficient 

Lift  coefficient  at 

Moment  coefficient 

Constant  defined  in  (9.5) 

Maximum  camber  of  the  hydrofoil 

Function  defined  in  (10.3) 

Integral  defined  in  (6.5) 

Integral  defined  in  (6.6) 

Constant  defined  in  (  6.4) 

Cavity  thickness  distribution  (see  Figure  1) 

Discrete  cavity  thickness  distribution  (see  also  C.l) 

Integral  defined  as  I/*j)  -  j ^  ^  ^ - -  oLi 

/  f  «  /-A  2* 

Integral  defined  in  (6.2) 

Integral  defined  in  (5.9) 

Number  of  uniform  intervals  to  apply  Simpson's  rule 
Integral  defined  in  (6.6) 
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€ 

Pco 

P(B) 

a 

r 

Sf^) 

i 

'It 

V 

v/x; 

v/x; 

V 

/^x,y; 

)^rx; 

z 


Cavity  length  (see  Figure  1) 

Free  stream  pressure  (see  Figure  1) 

Integral  defined  in  (A. 21) 

Cavity  source  distribution  on  the  cavitating  hydrofoil 
(see  Figure  1) 

:  Defined  as  ^ 2^  ) 

:  Quantity  depending  linearly  on  the  solution  (T^  ^ 

(i.e.  ^  j  i ) 

Constant  defined  as  /"=  ^  or  ^  {j.*  i’-y* 

Integral  defined  in  (6.7) 

Integral  defined  in  (A. 20) 

Constant  defined  as  i  = 

Maximum  thickness  of  the  hydrofoil 

X-  component  of  the  perturbation  velocity  (see  Figure  1) 
Free  stream  velocity  (see  Figure  1) 

y-  component  of  the  perturbation  velocity  (see  Figure  1) 
Induced  velocity  in  wake  by  the  vorticity  distribution 
Defined  as  V^X.)  =  V(x)  ) 

Cavity  volume  (including  the  foil) 

Coordinate  system  for  the  hydrofoil  (see  Figure  1) 

Camber  distribution  of  the  hydrofoil  (see  Figure  1) 

Thickness  distribution  of  the  hydrofoil  (see  Figure  2) 

Transformed  X  defined  as  ^rl/-— —  .*  0  ^  ^  "t 

r  X  ' 


Angie  of  attack 
Ideal  angle  of  attack 

Vortex  distribution  on  the  cavitating  hydrofoil 


rM 

ex 

e 


I 


(see  Figure  1) 


Defined  as  ) 

Increment  in  the  angle  of  attack  o< 


Similarity  parameter 


Tranformed  ^  defined  as  ^  /  ^■^’7  ***  ^ 

€r 

Distribution  of  the  lower  surface  of  the  hydrofoil 
(see  Figure  1) 

Transformed  defined  by  the  relationship: 

Defined  as  -  -i.  . 

<r  Bx 

Same  as  x 

Radius  of  curvature  at  the  leading  edge  of  the  hydrofoil 
Cavitation  number;  0's:  ) /fz 

vapor  pressure 
p>  :  density  of  the  fluid 
Integral  defined  in  (6.4) 

Transformed  ^  defined  by  the  relationship: 

y  O^ep^TL 

Same  as  » 


defined  as 


c  :  Refers  to  camber  problem  (see  Figure  3) 

if  :  Refers  to  thickness  problem  (see  Figure  3) 

0(  :  Refers  to  flat  plate  problem  at  angle  of  attack 

(see  Figure  3) 


(4.9) 


-  ^  .  f  — -  . /<at  -f-  ^>-7 ).  4  ^ )’ 

n.fr^l)  Jii-V  '  /  /  Px  /  / 


Equation  (4.9)  gives  in  terms  of  an  integral  of  the  slope 

dx 

of  the  lover  surface  of  the  hydrofoil  and  functions  of  .  The 
numerical  implementation  of  (4.9)  for  a  general  shape  hydrofoil  is 
described  in  §  10  .  Equivalent  expressions  for  but  also  for  ^  , 
and  have  first  been  given  by  Hanaoka  in  [  1  ]. 

In  this  work,  ve  focus  more  on  the  cavity  shape  than  on  the 
hydrodynamic  characteristics  of  the  cavitating  hydrofoil.  Thus  ve 
proceed  by  giving  general  expressions  for  the  cavity  volume \^and  the 
cavity  source  distribution  , 


(4.10) 


5.  CAVITT  VOLDME  V 


The  cavity  volume  V  vill  be  given  by  the  expression 

f  ^(x)’Ct-K  ST  X  •  ^^x)l  -  d£t.  .  clx  = 

^  ^  0 

^  /  X- .  Ux  =  -<r.  f  (5.1) 

o  ’  2^0  O'* 

vhere  ve  have  made  use  of  (1.4),  (1.10),  (3.1) 

Nov  plugging  (3.6)  in  (5.1),  by  changing  the  order  of  integration  and 
by  using  A. 11,  A. 12  and  (4.4),  ve  get: 


or  - 

2  ^ 


r 

°  . 


(4.5) 


using  the  transfomstion  (3.8). 

By  plugging  (3.10)  in  (4.5)  and  changing  the  order  of  integration 
ve  get: 

“  ^  /"  ;r  l/IEU.  [21  _  ,9^7  d.y} 

^  J  Lz 


‘-f  - 

--t-  -/w-f?  -  <>•]■  ih) . 

■t 


ft  2^  £-»7 


Where: 


ih)= f  - i_ .  ar? 

1  2r 


But 


••  rrv)  =  - 


n. 


f t- 


(see  A.  17) 


.  = 


(4.6) 


««■=  «=/p^vr^  ,4.71 


Therefore  the  closure  condition  (4.6)  becones: 

.t 

L 


(4.8) 


And  by  using  equations  A. 5,  A. 6  and  3.1  ve  finally  get: 
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Katta-condition  (3.4). 


By  plugging  (3.8)  in  (3.6)  ue  get: 


(3.11) 


The  ezpreeeionc  (3.10)  end  (3.11)  give  the  fonel  eolation  to  our 
problem  ee  soon  ae  ve  detcmine  the  correeponding  cavitation  noaber 


4.  CAVITATIOM  ITOMBER 


The  closure  condition  (3.5)  by  using  the  expression  of  j(x)tx<m 


(3.6),  becomes: 

f _ ^  • 

But:  flflElJ.  cCx.  =s  (see  A. 10) 

/  f  €-X  Z 


O  (4.1) 


Therefore: 


o 


(4.2) 


-  J-. 

Z  ^ 


~jdf  .  fj ^ 


by  changing  the^order  of  integration  in  the  double  integral  of  (4.2). 

But:  S =s /T  (see  A. 13a) 

JVZ-X  x-f 

c 


(4.3) 


(4.4) 


Z  z.fl  J 


;  o<x<l 


(3.3) 


(3.4) 


/ 


<9^  X  r:  o 


(3.5) 


Using  sinilsr  procedures  ss  in  Ix]  and  [2]  we  derive  the  solution 
to  the  system  of  singular  integral  equations  (3.2) -(3.5)  as  follows: 


Inverting  first  (3.2)  we  get: 


I  X'S 


notice  that  this  is  the  unique  solution  which  behaves  like  — at 
the  trailing  edge  of  the  cavity  (T.T.  Vu's  singularity). 


Plugging  (3.6)  in  (3.3)  we  finally  get: 

I 

O  ^ 


Using  the  transformation: 


e-y. 

(3.7)  becomes: 

L  r*  if, 


/  V=  — 


2n  J  fi’f- •>!*■) -(^->1)  4 

Inversion  of  the  (3.9)  gives:  . 


^  e-i 

t _ 

'■z‘-  Z  '  /4-Z*- 


yfi),  _  ^.(i* 


_ 


(3.6) 


(3.7) 


(3.8) 


(3.9) 


(3.10) 


notice  that  (3.10)  gives  the  unique  solution  which  satisfies  also  the 


'_U  •  -  •  •  •  -  m\  -  ta  •  I.  •  _  •  V  •  *  .  •  -V  -  •  JV  liV  _%  -  *  -•»  -N  _*»  _N  k_V  * 
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cavity  Icnctb  (aee  Fig.  3a) 


b)  Affine  ca»ber  ■caalinca  and  tbickneaa  foraa 


(2.8) 


where  Q.  and  Qt  correapond  to  easber/cbord  ratios 
r 

^  and  ^  respectively  and  aero  angle  of  attack 
(aee  Fig  3b). 


A  aiailar  relationship  is  also  true  for  affine  thickness  foras. 

Therefore  for  a  given  caaber  or  thickness  series  we  need  to  solve 

the  problcB  for  one  value  of  the  siailarity  paraaeter  and  one  angle  of 

attack. 

The  decoaposition  described  above  can  also  be  done  when  the  cavity 
does  not  start  at  the  leading  edge  of  the  foil.  In  that  case  both  the 
leading  edge  and  the  trailing  edge  of  the  cavity  have  to  reaain  fixed 
for  each  eleaentsry  problea. 


3.  gOLPTIOB  TO  THE  CEHEBAL  8UPERCAVITATIWG  PROBLEM 

It  is  convenient  to  non-diaenaionalise  equations  (1.6) -(1.9)  by 
dividing  by  C^'^)  and  calling 

I  ~  •ir' '  (3.1) 


)  ^,ir  / 


(T 


;  o<x<e 

2  2-r^J  /-X  2  J 


(3.2) 


with  ^ "" 


for 


^  O'  «S’  € 
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a)  Ca«ber  problcn  vith 

b)  Thickaets  problaa  vitb  ^ ^ — 

c)  Flat  Plata  at  angle  of  attack  o(. 


The  solution  vill  be  tbe  linear  superposition  of  the 


solutions  of  those  problems: 

(2.3) 

y(x')  =  riCx)  -h  Yi-M  -f-  Cx) 

(2.4) 

=r  %(x)  -h  ^c<(x) 

(2.5) 

Furthermore,  for  each  quantity  Q  which  depends 

solution,  we  will  have: 

linearly  on  the 

Cl  ^  +  Qo^ 

(2.6) 

For  ezaaple  can  be:  ^ 

Cl  ,  ^ 

SoBe  particular  cases  can  then  be  examined: 


a)  Change  of  anale  of  attack 

a  -  a'  -h  (2.7) 

where  Q  corresponds  at  angle  of  attack  o<  , 

In  other  words  if  we  change  the  angle  of  attack  by  A»<  then  the 
solution  for  tbe  same  cavity  length  will  change  by  an  amount  which  is 
the  solution  to  a  flat  plate  at  angle  of  attackid^  and  tbe  same 
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Figure  2 
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••  "-r 


V  VJ.' 


The  cavity  thickneas  ^0^)  and  the  cavity  caaber 
will  then  be  given  by  the  following  relationabipa: 


(1.10) 


CM  =  j 

z 


o 


(1.11a) 


?X 


d.  ^ 


V(x) 


{ 

2n.  J 


(l.llb) 


5  ^  <x<  C 


(1.12) 


y/(x)  i  induced  velocity  in  the  wake  by  the  vorticity  distribution  irfx) 


2.  LIHEAR  DECOMPOSITlOll 


If  the  hydrofoil  has  a  canber  aeanline  y^^x)  •  a  tbicknesa 
distribution  and  an  angle  of  attack  ,  then  the  lower 

surface  will  be:  (see  Fig*  2) 


(2.1) 


£5  =  ^  __  ^ 

<?x 


(2.2) 


Fron  equations  (1.6)-(1.9),  it  is  easily  seen  that  as  long  as  we 
keep  ^  fined  then  the  supercavitating  general  shape  hydrofoil  probli 
can  be  considered  as  the  superposition  of  the  following  three 
elcncntary  problems:  (see  Fig.  2) 


k(^) 


X  =  d. 


)  finite  at 


d)  Closure  condition 


Cavity  thickness  at  the  trailing  edge  of  the 
cavity,  ‘LfC')  -  O 


(1.4) 


e)  Condition  at  infinity 


sod  v-*  o  at  infinity 


(1.5) 


The  boundary  conditions  in  terms  of  the  unknown  distributions 
and  •  will  become: 


with  =  O  for  i  ^ 

^  z  /  r(i)-  _  ir  .  ^ 

Z  Z.n.  }  ^-x  ?x  ^ 


crx  <  I  (1-7) 


yfi)^  (p 


(1.8) 


j  Cj  (x")  X  - 


(1.9) 


The  boundary  condition  (1.5)  is  automatically  satisfied  by  the  use  of 
vortices  and  sources  for  the  representation  of  the  velocity  field.  For 
given  €  and  »  the  equations  (1.6)  -  (1.9)  give  us  the 

unique  solution  and  cj  , 


w  *\>  •> 


1.  FORMPIATIOW  OF  THE  PROBLEM 


ConsidcT  a  attpereavitating  hydrofoil  of  chord  1.  in  onifon  flow 
“IT  and  aahient  preaaure  p  (aee  Fig.  1.).  Given  the  cavity  length  ^ 
and  the  lover  aurface  the  hydrofoil  we  want  to  determine  the 

corresponding  cavitation  number  ff"  and  the  shape  of  the  cavity. 

According  to  the  linear  theory  the  perturbation  velocity  field 
(components  u  and  v)  due  to  the  cavity  and  the  foil  can  be  generated  by 
a  distribution  of  vortices  ^(  >c  )  and  sources  (X  )  »  along  the  slit 
from  X~0  to  (see  Fig.  1.). 

The  linearised  boundary  conditions  are  (throughout  our  analysis, 
ve  assume  that  the  cavity  starts  at  the  leading  edge  of  the  foil): 


a)  Dynamic  boundary  conditions 


/  o<x<  €  ^ 

(1.2 

;  icx<-e  , 

II 

1 

(1.: 

h)  Kinematic  boundary  condition 

II 

(1.: 

c)  Kutta  condition 


compact  than  the  ones  given  by  the  previous  authors.  Be  also  gave 
series  representations  for  his  results  when  the  hydrofoil  shape  could 
be  expressed  in  terms  of  polynomials  of  the  chordvise  coordinate. 

In  this  work,  the  linear  problem  is  stated  in  a  simple  way  but 
equivalent  to  the  one  described  in  [1],  in  terms  of  integral  equations 
of  cavity  source  and  vortex  distributions.  The  principle  of 
decomposition,  first  introduced  by  Fabula  for  a  hydrofoil  with  cut-off 
ventilated  trailing  edge  [111  and  for  a  super cavit at ing  hydrofoil  at 
different  angles  of  attack  [12],  is  extended  to  our  problem  which  is 
divided  into  the  camber,  the  thickness  and  the  angle  of  attack 
problems.  Formulas  for  the  cavitation  number,  the  cavity  volume  and  the 
cavity  source  distribution,  are  given  in  terms  of  integrals  of  the 
lower  surface  of  the  hydrofoil  weighted  by  known  functions.  These 
integrals  are  computed  numerically  using  Simpson's  rule  with  respect  to 
an  appropriate  transformation  variable,  thus  avoiding  the  involved 
singularities  of  the  integrands.  The  cavity  shape  is  also  computed 
numerically.  The  results  for  a  flat  plate  and  an  elliptic  foil,  for 
which  analytical  solutions  are  available,  seem  to  be  very  accurate 
even  for  a  small  number  of  the  used  intervals  in  the  numerical 
integrations.  The  cavity  shape  for  a  flat  plate  is  compared  to  the 
nonlinear  result  given  by  Uhlman  [5],  and  is  found  to  be  very  close 
except  at  the  trailing  edge  of  the  cavity. 

Finally,  the  proposed  method  is  applied  for  the  KACA  a  "  0.8 
meanline  series  and  the  BACA  Four-digit  Wing  Sections  thickness  forms 
for  which  tables  and  some  cavity  plots  are  given. 


IHTRODUCTIOM 


The  problem  of  a  eupercavitat ing  hydrofoil  in  uniform  flow  has 
received  a  lot  of  attention  in  the  past  three  decades,  individually  as 
well  as  a  first  step  towards  the  understanding  of  the  pbenoinenon  of 
cavitating  marine  propellers. 

The  non-linear  formulation  of  the  problem  is  not  unique  due  to 
the  variety  of  the  cavity  termination  models.  Furthermore,  the 
non-linear  problem  is  very  difficult  to  deal  with  analytically, 
especially  for  general  shape  hydrofoils. 

Linear  theory  was  first  applied  by  Tulin  [8]  to  the  problem  of 
a  supercavitating  flat  plate  at  small  angles  and  arbitrary  cavitation 
numbers.  It  was  subsequently  extended  by  T.T.  Wu  [9],  Geurst  [3], 
Parkin  [10]  and  Fabula  [12]  for  a  supercavitating  hydrofoil  of  general 
shape.  These  authors  have  worked  by  means  of  the  complex  velocity 
function,  and  have  given  expressions  for  the  cavitation  number,  the 
hydrodynamic  coefficients  and  the  cavity  volume,  in  terms  of  integrals 
of  known  quantities. 

Hanaoka  [1]  stated  the  linear  problem  in  terms  of  integral 
equations  of  unknown  distributions  of  sources  and  dipoles  on  the  foil. 
Inversion  of  these  integral  equations  produced  integral  representations 
for  the  cavitation  number,  the  hydrodynamic  coefficients  and  the  slope 
of  the  upper  and  lower  cavity  surfaces  in  terms  of  the 
shape  of  the  hydrofoil.  His  formulas  are  less  complicated  and  more 


Casting  (5.2)  in  terns  of  (3.8)  we  have: 


(5.8) 


where 


'fi- 

Now  ^(^'j  c&n  be  expressed  in  terms  of  y  and  as  follows 


(5.9) 


(see  A. 18) 


! 


%f6  I  J 

•f  V  t (5.10) 

And  finally  the  cavity  volume  is  given  in  terns  of  is.  .  dr 
and  «  as: 

if[v-  »■  •'7/ 


(5.11) 
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Starting  fron  equation  (3.11)  and  plugging  in  the  expression  for  JT 
from  (3.10),  we  get:  , 


-22- 


(6. 


=-  -Ir- 


cj) 


We 


call  the  double  integral  in  (6.1)  • 


-cj 


0 


.CO 


where: 


(i-t^^) 

Hov,  to  change  the  order  of  integration  in  (6.3)  ve  apply 
Poincare-Bertrand'a  formula: 

.-b  A 


J(i)=  ( - -dv,  •  if,  ir 


l(i)=  \<iu  :  ici 


We  then  focua  on  the  integral 


l<[i)^  f  dco  i  - .tbyj  r- 

and  ve  call :  ^  ^  1  /  j 


(6.2) 


(6.3) 


(6.4) 


(6.5a) 


(6.5b) 


(6.6) 


(6.7) 


O’  ^  •  f 

But/?/^^  can  be  cxprcaaed  explicitly  aa:  (aee  A. 19) 


Hfu,  2;  =  - 


7L 


?;  =  — 1L 


/t 


f¥-) 


% 


; 


(6.8a) 

(6.8b 


Plugging  (6.5)  in  (6.1)«nd  uciog  (6.6),  (6.8),  (A. 20)  and 
(A. 21)  va  get: 


~±- .  ff^T- 


J— 

Z4UJ 


^  <  -t 


(6.9a) 


71  '  '  /  ?■  /  Jf 


— i - / /rVi  1 2  -h  /r^i  ^  j  - 


<  2-  /  4?.i£!ri 


-  f  Y" -w;  V/^- sfe  •  ^  ^ 


for 


(6.9b) 


Integrating  (1.10)  and  naing  (3.1)»  ve  get: 

c 

For  the  noaerical  integration  of  (7.1)  aee  §1.0. 


(7.1) 


8.  CAVITY  CAMBER  DISTRIBUTIOM  COt> 


The  camber  meanline  of  the  cavity  and  foil,  d:  (x),  for 
it  explicitly  determined  from  (1.11a)  in  terma  of  the  foil 
characteriatica  and  the  cavity  thickneaa  ¥otI<^X<^  it  ia 

going  to  he  determined  from  (1.11b)  after  we  firat  ezpreaa  V(x)  in  terms 

of  ia  and  ^  : 

?X 

From  the  definition  (1.12)  ve  get: 


(8.1) 

(8.2) 


Expreaaing  (8.1)  in  tema  of  the  trana  format  ion  (3.8)  and  then  naing 


(8.3) 


Using  now  A.22»  A. 20.  A. 21,  ve  get: 


_i_.//+2*o.  /  JIKT .  . -J— .  du) 

Js  t.n.  V  J  j'j-t-u  g.co 

0 

for 


end  by  integrating  (1.11b)  ve  get: 


CCx)=:  C(l)  +  <r.  rY(g).dS 


with 


cfi)  =  f /i;  ^ 


(8.5) 


(8.6) 


(8.7) 


9.  AHALYTICAL  80LUTI0WS  FOB  SOME  HYDKOFOTT.g 


9.1  Flat  Plate  (see  Fig.  4a) 
For  this  case  ve  have 


where  o(  ■  angle  of  attack 


26- 


'A  pj, 


And  tbe  •zprestion*  (4.9),  (6.9),  (8.5)  by  vting  tbc  fonmlas  in 
Appendix  A  reduce  to:  (tee  also  [  ^  1  and  [  3  ]) 

^=:  ^  <9-2) 

=  fc- y  Z<-t  {9.3a) 


(9.3b) 


(9-4) 


(9.5) 


Also  by  using  (5.11)  and  (9.1)  ve  Bust  recover  tbe  result  by  Ceurst  in 
I  3  ]  for  tbe  cavity  voluae  : 


9.2  Elliptic  Foil 

Consider  an  ellipse  of  seai-azes  )  (>••  fig*  4b)  at  unifons 

inflow  Z/  .  According  to  tbe  linear  theory  tbe  perturbation  velocity 
on  tbe  ellipse  will  be  constant:  (see  (  4  ]  ) 

(9-7) 

Tberefore,  if  ve  consider  tbe  lover  left  part  of  tbe  ellipse  as  a 
super cavitating  hydrofoil  at  cavity  length  2 X chord  then  tbe  shape  of 

-28- 


.N  .  V 


I.J  1,1.  J...  II lH..i.iV-l«'i 


the  tmyity  will  be  the  rett  of  tbe  ollipoe  (in  liooor  tboery  olvajs)* 
and  w€  will  alao  get  by  iMpeetion: 


o’s  2.  e 

(9.8) 

y“=  /!•  £ 

(9.9) 

-  'll  > 

(9.10) 

••N 

ii 

lit 

2^  ir 

(9.11) 

Bqoatioas  (9.8)-(9.11)  cao  alao  be  recovered  by  plogBin^  in  (4.9), 
(5.11),  (6.9)  and  (8.5)  the  expreaaion  for  : 


•  0<€O<i 


(9.12) 


10.  BUHERICAI  AHALTSIS 


For  general  abape  bydrofoila  ve  compote  tbe  derived  fomolaa  in 
tbe  previooa  aectiona  nnnerieally. 

Vbenever  an  integration  in  ^  from  O  to  ia  involved  ve  amke  tbe 
tr ana format  ion: 

Slyi^l ^/Z')  (10.1) 

and  then  oe  nae  fimpaon'a  mle  vitb  K  uniform  intervals  li.e. ,  2-K+ 1 

pointa). 

For  asample  equation  (4.9)  becomea: 


ITS 


■h  L  t.  si  >i*-(9/2.)j 
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(10.2) 


The  integrand. 


(10.3) 


bebavea  at  the  liaita  of  tbe  integration  as  follows : 
at  ll^s  O  : 


If  the  foil  has  leading  edge  radiua  f  then: 

fj ~  '/^vf •/  '  »s  /->° 

and  _  _ 

px  U  fT 

or  by  using  (3.8)  and  (10.1) 


(10.4) 


_ .rn  i_,j_  ^ 

i2.e  ^  *  j4 


(10.5) 


and 


(10.6) 

at  i^sA*  linite  as  long  as  is  (which  is  the  usual 

case). 

Therefore  by  working  in  tens  of  we  avoid  tbe  square  root 
singularities  in  (4.9). 


We  apply  the  sane  technique  for  the  integrals  involved  in  the 
fonnlas  (5.11),  (6.9)  and  (8.5). 

To  find  the  shape  of  the  cavity  we  nuaerically  coapute  the 
fonnlas  (7.1)  and  (8.6)  by  using  Siapson's  rule  in  the  transfoned 

variable  ^  :  (see  also  Appendix  C) 

/=  (10.7) 


<r  6  •< - >  o  ^  4>  ^  n. 
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To  coapate  (7.1)  ve  need  to  know  hov  behaves  at  the  leading  edge 
and  trailing  edge  (?c«e)  of  the  cavity.  Using  fomulas  (6.9)  ve 
get  the  following  asyapotic  behaviors:  (see  Appendix  U) 
at  : 

which  states  that  the  cavity  for  roonded  nosed  hydrofoils  starts 
tangent  to  the  upper  part  of  the  foil. 


(10.9) 


which  is  the  expected  square  root  singularity. 

The  integral  in  (10.9)  is  coaputed  nuaerically  after  asking  the 
transforaation  (10.1). 

The  sensitivity  analysis  of  the  involved  nuaerical  integrations 
for  the  special  cases  of  a  flat  plate  and  an  elliptic  foil  is 
described  in  Appendix  C. 

A  coaparison  for  the  cavity  shape  of  the  linear  to  the  non-linear 
theory  as  developed  in  [  5  ]  can  be  seen  in  Fig  5  .  Ve  observe  that 
the  two  shapes  are  very  close,  except  at  the  trailing  edge  of  the 
cavity,  which  is  doe  to  the  iaposed  tiabooshinski  aodel  in  the 
non-linear  theory. 

Finally,  the  analysis  described  above  is  applied  to  the  series  of 
NACA  a  -  0.8  meanlines  and  NACA  00  thickness  forms  for  which  tabulated 
results  along  with  some  cavity  plots  are  included  in  Appendix  D. 
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APPENDIX  A 

LIST  OF  INTEGRALS 

In  this  Appendix,  a  list  of  integrals  used  throughout  our  analysis 
is  given,  along  with  instructions  for  their  derivation. 

The  first  integrals  (A.l)  ■*  (A. 9)  have  been  taken  from  Appendix  A  of 
[  2  1 
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=  -  o.s.fo.  o.iz^oo.x  -  0.5S16  0  -f 


4-  (?.Zff^30.  X  -  o.4oi 


S o  .  ) 

'  J 


for  O  ^  X  <  1 


The  results  for  Of  and  T^are  given  in  Tables  9  and  10 
(see  also  Figures  10  and  11)  along  with  some  cavity  plots  in  Figures 
7a,  7b,  7c,  8a,  8b,  8c.  The  values  for  ^  and  l^jhave  been  computed 
with  =  50  in  order  to  get  the  first  five  decimal  places  correct 
(although  we  get  the  same  accuracy  with  much  smaller  for  larger 
and  the  cavity  plots  have  been  made  with  K  =  30  points  (although  the 
cavity  shape  has  been  found  very  insensitive  to  for  a  broad 
range  of  ^  ).  Most  of  the  plots  have  no  physical  meaning  but  if 
combined  with  the  appropriate  angle  of  attack  so  that  the  cavity 
clears  off  the  upper  part  of  the  hydrofoil  then  they  give 
meaningful  solutions  to  our  problem. 

In  case  our  hydrofoil  consists  of  a  HACA  a  0.8  meanline  of 
maximum  camber  ^  ,  a  NACA  00  thickness  form  of  maximum  thickness 

operates  at  an  angle  of  attack  a(,  ,  then  according  to 
Section  2,  any  quantity  Q,  can  be  expressed  as:  (keeping  always  the 


same  cavity  length) 


L  / 


Q-~ 


A  •  ^  _  COlO  ^ 


Aof^j.c  aio 


r 


Q 

r 


where: 


(^o.e 

^90lO 


:  Q,  for  HACA  a  -  0.8,^^  -  0. 


0679 


(D.4) 
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SUPERCAVITATING  NACA  a  »  0.8  MEANLINES  AND  NACA  00  THICXNESS  FORMS 


The  analysis  described  in  the  previous  sections  has  been 
carried  out  for  a  NACA  a  *  0.8  meanline  (with  “  0.0679)  and 

NACA  0010  thickness  form. 

For  the  NACA  a  *  0.8  at  ;  we  put  as  the  analytical 

expression  from  [7]: 


.  1 .  (I- x;  ^  X.  x/-  X  /■«  -  x/y  - 

for  0  <  X  <  i 


/ 


where : 

a  •=  C.8,  C.  =  1,  *  1.5A* 

i  * 

I  -  -  0.009297,  ^  =  -  0.3039 


For  the  NACA  0010  thickness  form  at  ^  0*  we  also  take  from  [7]: 


(D.l) 


(D.2) 


-48- 


ELLIPTIC  FOIL  AT 


e  =  2  (  £  =  1) 


Table  7 


CAVITY  THICKNESSES  A-i 

's  (  k:  =  20) 
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Numerical 

Analytical 
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^  =  2  (£  =  1) 


Table  6 


CAVITY  THICKNESSES  (  K  =  5) 


. . 

Nvunerical 

Analytical 

r,  ^  r-  ; )  f- 
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CAVITY  THICKNESSES  A;  *S  I  K  =  10) 
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000000 
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CAVITY  SOURCE  DISTRIBUTION  ^ {X. 


Numerical 


Analytical 


2.06474 


1.33333 


0.75000 


0.57735 
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0.00000 


k  =  10 


2.06474 


1.33333 


0.75000 


2.06206 


1.33270 


0.74987 


0.57735 

0.57727 

0.31449 

0.31444 

0.10050 

0.10046 

0.00000 

-  0.00005 

FIAT  PIATE  AT  ^/c  =  1.4 »<  =  1  rad 


Table  3 


CAVITY  SOURCE  DISTRIBUTION 

Vc 

Analytical 

Numerical 

0. lOOOOOOOD-02 

0.289485380+01 

0,289485380+01 

■ 

0* lOOOOOOOD-Oi 

0,160065280+01 

0.160065280+01 

0.  lOOOOOOOEifOO 

0,830929700+00 

0,830929690+00 

0,3000000011+00 

0,541655360+00 

0,541655360+00 

0,5000000011  +  00 

0,394577540+00 

0,394577540+00 

0,800000000+00 

0,212297250+00 

0,212297250+00 

0,100000000+01 

0,692101400-01 

0,692101400-01 

0, 120000000 » 01 

-0,699557530+00 

-0,699557530+00 

0,135000000+01 

-0,202597900+01 

-0,202597900+01 

0,139900000+01 

-0,155899890+02 

-0,155899890+02 

FLAT  PLATE  AT  «  =  1  rad 


Table  1 


Analytical 


CAVITATION  NUMBER  g* 

Nvunerical 


1C  =  20 


1C  =  10 


JC  =  5 


1.1 

6.324555 

6.324555 

6.324560 

6.329471 

1.3 

3.651484 

3.651484 

3.651484 

3.651730 

1.5 

2.828427 

2.828427 

2.828427 

2.828405 

1.8 

2.236068 

2.236068 

2.236068 

2.236057 

2.0 

2. 

2.000000 

2.000000 

1.999995 

Table  2 


CAVITY  VOLUME 


4 

Numerical 

/c 

Andly^xcdX 

1C  =  20 

H 

O 

k  =  5 

1.1 

1.272656 

1.272656 

1.272668 

1.282604 

1.3 

1.327715 

1.327715 

1.327715 

1.328218 

1.5 

1.554475 

1.554475 

1.554475 

1.554295 

1.8 

1.975729 

1.975729 

1.975729 

1.975635 

2.0 

2.288818 

2.288818 

2.288818 

2.288768 

’my. 


For  a  aupercavitating  flat  plate  at  K  -  4  and  different 
cavity  lengths,  the  values  for  ^  and  H^are  given  in  Table  8  (see 
also  Figures  10  and  11),  and  some  cavity  plots  are  shown  in  Figures  6a, 
6b,  6c. 


•  •  •  • 

•  •  •  a  ■  ' 


The  numerical  analysis  as  described  in  Section  10  has  been 
tried  for  the  special  cases  of  a  flat  plate  (Tables  1  and  2)  and  an 


elliptic  foil  (Tables  4-7)  for  which  the  analytical  results  are  given 
in  Section  9.  The  results  are  very  good  even  for  small  . 

For  the  computation  of  the  cavity  source  distribution  frx>  for  a  flat 
plate  *  50  has  been  used  for  0.<'>'<'0.1  and  1.<'>«^1,2,  and  J<  »  10 
elsewhere,,  to  give  us  the  excellent  results  as  shorn  in  Table  3. 

The  cavity  thicknesses  have  been  computed  using  (7.1)  and  the 


transformation  (10.7)  as  follows; 


=  *-i 


cL^ 


■%  eC^  2  ^ 


(C.l) 


(C.2) 

(C.3) 


The  integral  in  C.l  has  been  evaluated  by  using  Simpson's  rule  in  one 
interval  between  and  d,  A  check  to  our  numerical  integrations 

•  I 


should  be  the  closure  condition; 


=  O 


(C.4) 


A  similar  (if  not  better)  convergence  is  found  for  the  camber 
distribution  C  (.X)  in  the  wake  of  the  cavity. 


-41- 


(B.IO) 


as  ^-^09  and  for  both^eitfaer  sharp  or  rounded  nosed  hydrofoils. 


-40- 


Ffi)  =  J—  .  ■  d-co 

a-i-eo 


(B.5) 


Sine.  no.  tbon  ^.0—'. 


as  C0~*0  and  P(i)  ^  as  Z-*0  ,  where; 

A 


f^(2h  - 


■4 


Z.G- 


fP  J  fP.(2ri-t*>)  ifp./F 


.  ^rc 


-•> 


Thus: 


F(2)  ^ 

(f./F 


2‘  O  .  Jk,  —  ^  /r.  ^ 


as  ^  ® 


Z  -*0 


(B.7) 


and  finally  — 


- 2.  ^7^;  as  2  ^  ^  (B.8) 

which  means  that  for  round  nosed  hydrofoils  the  cavity  has  the  same 
slope  and  curvature  with  the  wetted  part  of  the  hydrofoil  at  the 
leading  edge. 


b)  Sharp  nosed  hydrofoils 


The  asymptotic  behavior  of  will  be;  (as  jj  — >  ^  ) 

since  ^{Z)«  ^  -  as  Z-^O  (even  for  NACA  a  -  0.8  meanline) 

Now  at  the  trailing  edge  of  the  cavity  (j?=*^  or  ) 

by  using  (6.9b)  we  get  the  asymptotic  behavior  of  ^fZ)  ass 


(B.9) 


-39- 


APPEKDIX  B 


ASYMPTOTIC  BEHAVIOR  OF  THE  LEAT^TMC  AMD  TRAILING  EDGES  OF 


THE  CAVITY 


First,  at  the  leading  edge  of  the  cavity  and  foil 
(  Z  “  0  or  X  =  0)  we  consider  the  following  cases: 
a)  Round  nosed  hydrofoils 

If  the  leading  edge  radius  is  jP  we  have 

^  -  )fL^.K  as  X  <9 

and  by  using  (3.1)  we  get: 

But  i  ^  as  ,  according  to  (3.8) 

rr 


(i)  ^  - 


To  find  the  asyiiq>totic  behavior  of  (6.9a)  as  2-i^O  we  need  first  to 
work  on  the  integral: 


WWWW! 


and  finally  «e  get: 

RftOjir)  =r - ^ - 

(i+to) 


2^i,  ^ 


'ii-f 


!  O 


(A. 19a) 


^)  =  -  -JL—  ./41L  — ^ — JiH:  ,. 

Sinilarly  as  in  the  previous  sections  ve  can  derive  : 

^  _ _ 

5/?;=  /  .  -1 - f£_ .  = 


(A. 19b) 


-L..  -5^  ./^_  _2_ 

L  J  f i-i-i 


•  2>0  (A.2C 


P(i)=  f  /_*L-'  - = 

J  h-‘^  1+u^  B-x) 


= //■^-  — 2 - = 


c 


1+2^ 


^-7 


where: 


A-  V  c  =  — 

/  «  -  "  77^  /  ^ 


•  ’  •  r 

1  finally  by  plugging  in  the  known  integrali  we  get: 

- - — -  •  f aCi- 


(A.l- 


where : 


cL  =  ^r^-hd.  '  /  ^ 

'l-gl-  M-i '  /  /,  _ 


=  / i/^z_? .  - j  O  ^-y]  ^  ir 

Jl  -^-2=  ^ 

Similarly,  we  break  the  integrand: 

_  A^2-hB  ^  CZ-hO  ^  e 


C-I  *>  _*■  j  ^  /7«  *  ^ 

i-i-z 


where: 


A=  JhinzIL  ,B  =  £tLl±L 

At*  7^/*^ 


r-  J±2 


and  finally: 


ty-il  E-  vh-il 


»  /  7  - 


7/^)=  ^ - rj.fqrl- 

-  -f- 

=)  M.fcOj  B)  =  flfHE! ,  ^ - L^.dri  £> 


Ve  have: 


Q  • 
V  * 


o<>  ={.S4  : 


^  for  flat  plate  at  angle  of  attack 
-  4® 

ideal  angle  of  attack  for 
HACA  a  »  0.8  =  0.0679) 


For  example.  If  our  hydrofoil  is  a  combination  of  a  RACA  a  ■  0.8 

meanline  with  4  *  0.0679  and  a  RACA  0010  thickness  distribution 

0 

and  operates  at  an  angle  of  attack  ^  1  then: 


n  /T*  -I.  O'  .  ^ ^ 

”  0.9  Oo£o  ^ - 


V'  V  V  ^  Y 

^  ^0010  ^  ^  "4* 


(D.5) 


(D.6) 


If  ^  =  1.4  then  by  substituting  in  (D.5)  and  (D.6)  the  appropriate 
values  from  Tables  8,  9  and  10  we  get: 

^  -  0.52946  and  Y  “  0.21782 


Also  for  the  cavity  thicknesses  we  will  have: 


(D.8] 


which  is  il''ustrated  in  Figure  9 
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950 

0.62443 


A 


0.36051 


0.27925 


,.492 


0 . 2360 1 


0.22077 


0.20814 


0.19746 


0.18827 


0.18026 


0.17319 


0.16689 


0.16123 


0.15611 


0.15145 


0.14718 


0.14325 


y/c* 


0.10030 


0.03885 


0.08672 


0.08753 


0.08972 


0.09269 


0.09618 


0.10003 


0.10417 


0.10852 


0.11307 


0.11778 


0.12263 


0.12761 


0.13272 


0.13793 


0.14325 


0.14867 


0.15418 


2.000 


0 . 139  63 


0.15979 


OOo 


26 


0.04355 
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